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Abstract 



We show that the ID Courant-Snyder theory can be con- 
sidered as a partial case of the multidimensional theory of 
coupled particle beams with equal eigenemittances. 

INTRODUCTION 

The parametrization of coupled beam motion has been 
studied intensively over the past decades. Nevertheless, 
there is still no representation of general coupled motion 
that would be as complete and as widely accepted as the ID 
Courant-Snyder theory. But is it really so, that the Courant- 
Snyder theory is, essentially, the theory which is applicable 
only to the ID linear motion? Our answer is that the one- 
dimensionality of motion is not the main source of the ele- 
gance of the Courant-Snyder approach. It is the property of 
the beam matrix to be proportional to the matrix which is 
simultaneously symmetric positive definite and symplectic. 
Because this property is independent from the phase space 
dimensions and is the characteristic property of the parti- 
cle beams with equal eigenemittances, the purpose of this 
paper is to show that the ID Courant-Snyder theory can 
be considered as a natural partial case of the multidimen- 
sional theory of coupled particle beams with equal eigene- 
mittances. It is also possible to say (in the opposite way), 
that in this paper we extend the ID Courant-Snyder for- 
malism to the multidimensional theory of beams with equal 
eigenemittances. Due to space limitation, we consider only 
the most basic theoretical questions and the more detailed 
study will be presented in a separate publication. 

BEAM MATRIX AND 
EIGENEMITTANCES 

Let us consider a collection of points in 2n-dimensional 
phase space (a particle beam) and let, for each particle, 

z = iq,p) T = (qi,---,qn,pi,---,Pn) T (l) 

be a vector of canonical coordinates and momenta. Then, 
as usual, the beam (covariance) matrix is defined as 

^ = ((z-(z)).(z-(z)) T )^(^ ^ Y (2) 
\ / y ^pq ^pp j 

where the brackets ( ■ ) denote an average over a distribu- 
tion of the particles in the beam and n x n submatrices of 
the 2n x 2n matrix £ satisfy 

^<?<? = ^<?<?' ^PP = ^P<? = ^gp" 

One says that the beam is uncoupled if all four submatrices 
of the matrix £ are diagonal matrices. By definition, the 
matrix £ is symmetric positive semidefmite and in the fol- 
lowing we will restrict our considerations to the situation 
when this matrix is nondegenerated and therefore positive 



definite. For simplification of notations and without loss of 
generality, we will also assume that the beam has vanishing 
first-order moments, i.e. (z) = 0. 

Let s be the independent variable (time or path length 
along design orbit) and let M(t) be the linear transfer ma- 
trix which propagates particle coordinates from the state 
s = to the state s = r, i.e let 

z(t) = M(t)z(0). (4) 
Then the beam matrix £ evolves between these two states 
according to the rule 

£(r) = Af(r)£(0)Af T (T). (5) 
If the matrix M(r) is symplectic and satisfies the relations 

M t (t)JM(t) = M{t)JM t {t) = J, (6) 

Whel ' e r ( I\ 

J ={-I J (7) 

is the 2n x 2n symplectic unit matrix and / is the n x 
n identity matrix, then the transport equation (0 can be 
transformed into the following equivalent form 



(£J)(t) = M(r) • (£J)(0) • M~ x (t). 



(8) 
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From this form of the equation (|5} we see that the eigenval- 
ues of the matrix £ J are invariants, because ((H) is a sim- 
ilarity transformation. The matrix EJ is nondegenerated 
and is similar to the skew symmetric matrix E 1 / 2 J E 1 / 2 

EJ = E^^E^JE 1 / 2 )^- 1 / 2 , (9) 
which means that its spectrum is of the form 

± iei, . . . , ±ie n , (10) 
where all e m > and i is the imaginary unit. 

The quantities e m are called eigenemittances and are 
generalizations of the usual ID rms emittances 

£m = \J {qL){p 2 m ) ~ (qmPm) 2 (11) 

to the fully coupled case JTJ. It is not difficult to prove that 
the set of all eigenemittances {e m } and the set of all pro- 
jected emittances {e m } of a given beam matrix E coincide 
if and only if the beam is uncoupled. 

The other approach to the concept of eigenemittances is 
the way pointed out by Williamson's theorem (see, for ex- 
ample, references in JT]). This theorem tells us that one 
can diagonalize any positive definite symmetric matrix E 
by congruence using a symplectic matrix T 

TET T = D, (12) 
and that the diagonal matrix D has the very simple form 

D = diag(A,A), A - diag(e 1; . . . , e„) > 0, (13) 
where the diagonal elements e m are the moduli of the 
eigenvalues of the matrix E J. The matrix T in ( TT2l is not 
unique, but the diagonal entries of the Williamson's normal 
form D (eigenemittances) are unique up to a reordering. 



TWISS PARAMETERS 
OF THE BEAM MATRIX WITH EQUAL 
EIGENEMITTANCES 

Let us assume that the beam matrix £ has all eigenemit- 
tances equal to each other and equal to the value e > 0. 
Then, according to the Williamson's theorem, there exists 
a symplectic matrix T such that 

S = eT- x T- T = eW. (14) 

The matrix W in ( TBI is independent from any particular 
choice of the diagonalizing matrix T in (TT~2t and is simul- 
taneously symmetric positive definite and symplectic. We 
will call it the Twiss matrix and will parametrize it as fol- 
lows 

/3 —a 

— a T 7 



W = 



(15) 



where the n x n submatrices (3 = (3 T , a and 7 = 7 T are 
the natural matrix generalizations of the corresponding ID 
scalar Twiss parameters. Due to symplecticity of the matrix 
W the matrix Twiss parameters satisfy the relations 



07 = I + c 

a f3 = /3 a 1 



7 a 



a T 7. 



(16) 
(17) 



(18) 



Because the matrix W is positive definite, both its subma- 
trices /3 and 7 are also positive definite and, therefore, non- 
degenerated. From this it follows that not all relations (|T6b - 
( |T8l are independent. For example, the relations (TTTj or the 
relations ( fT8l can be omitted. 

The inverse of the Twiss matrix W is given by the for- 
mula 

n-' = ( ; * ) . .19) 



7 

Q 



and, therefore, the natural multidimensional analogy of the 
positive definite ID Courant-Snyder invariant 

-1 _ 



(20) 



can be written as follows 

I cs = q T "fq + q T a T p + p T aq + p T /3p. (21) 

In general, the eigenemittances of the matrix £ have to 
be found as solution of the eigenvalue problem for the ma- 
trix X J. But if one wants only to know if all of them are 
equal to each other or not, then the solution of the eigen- 
value problem is not necessary. It can be done by simple 
matrix multiplication as explained in the following propo- 
sition. 

Proposition 1 The beam matrix £ has all eigenemittances 
equal to each other and equal to the value e > if and only 
if the equality 

(£ J) 2 + e 2 / 







(22) 



holds, i.e. if and only if the matrix (£ J) 2 is a negative 
scalar matrix. 



PARAMETRIZATION OF BEAM 
TRANSFER MATRIX, NORMALIZED 
VARIABLES AND PHASE ADVANCES 

Due to nonuniques of the diagonalizing matrix T in the 
Williamson's theorem the relation (fl4l can be considered 
as a multi-valued function which maps any particular ma- 
trix £ (and/or W) into some subset of symplectic matrices. 
Let us select some single-valued branch of this function 
which associates one, and only one, output to any particu- 
lar input. We will call any such branch as T-algorithm and 
the examples of the T-algorithms will be given below. So, 
let us fix some particular T-algorithm and let us substitute 
the representation (fl4l into the equation (0. Then, after 
some straightforward manipulations, we obtain 

(T(t)M(t)T -1 (0)) ■ (T(t)A/(t)T" 1 (0)) T = I, (23) 

which means that the 2n x 2n matrix 

R(t) = T(t)M(t) T -1 (0) (24) 

is orthosymplectic (i.e. orthogonal and symplectic simulta- 
neously). The equality d24l i. when written in the form 

M(r) = T-\t)R(t)T(0), (25) 

gives us a (familiar in ID) parametrization of the beam 
transfer matrix M(r), and if we will introduce normalized 
variables z n by the equation 

z(s) = T- 1 (s)z„(s), (26) 

then the dynamics in the normalized variables 

* n (r) - R{r)z n {Q) (27) 

is simply a rotation and the multidimensional Courant- 
Snyder invariant (fJT) takes on the form I cs = zj ■ z n . 

Although the motion in the normalized variables (l27l is 
not, in general, uncoupled motion, but from the point of 
view of the theory of beams with equal eigenemittances no 
additional simplifications are required, because 

En = (zn-zl) = el (28) 

is already a diagonal matrix. 

Let us turn now our attention to the concept of phase ad- 
vances, which does not have an unique choice even in the 
ID case (see, for example, discussion in J2])- The phase 
advances are the quantities which should be associated with 
the eigenvalues of the matrix R. This matrix is orthosym- 
plectic and can be partitioned into the form 

where the n x n submatrices C and S satisfy 

CS T = SC T , CC T + SS T = I. (30) 

All eigenvalues of the matrix R(t) lie on the unit circle in 
the complex plane, i.e. are of the form 

exp(±i/zi(r)), . . . , exp(±i/x n (r)), (31) 

and /i m (r) are the quantities which we will interpret as 
(fractional part of) phase advances. 

If the beam transport in (0 is periodic (i.e. if E(r) = 
E(0) and, therefore, T(r) = T(0)), then the equality (|25) 



tells us that the eigenvalues of the matrix R(t) are the same 
as the eigenvalues of the matrix M(r). It means that for the 
periodic beam transport the phase advances are uniquely 
defined independently from any particular choice of the T- 
algorithm. It is a very pleasant property, but it seems that 
it is the only property of the phase advances which does 
not depend from the choice of the T-algorithm. Let us, for 
illustration, consider three ID T-algorithms defined by the 
requirement for the matrix T in ( TBI to be in one of the 
following special forms 



T 



Ti 



T 3 = 



a/VP VP 



T 9 



^03 + 1) + ( 7 + l) 



Vy a IV^ 

1/V7 



7 + 1 a 
a P + l 



, (32) 



(33) 



Note that here T\ is the original Courant- Snyder choice and 
T3 is symmetric (and symplectic) positive definite square 
root of the matrix W~ Y . In all these cases the matrix 
R(t) = R(fj, m (r)) in d25T l is given by 



R(Hm(T)) 



cos(^ m (r)) sin(jU m (r)) 



(34) 



-sin(/i m (r)) cos(/U m (r)) 

and, in order to see more clearly the difference between the 
behavior of the phase advances p, m , let us assume that the 
dynamics is derived from the Hamiltonian 

H(r) = (1/2). (pl + k(r)ql). 

Then the phase advances \i m obey the equations 

dfii 1 d[ii k 
dr P dr 7 ' 



(35) 



(36) 



dp* 
dr 



k(P + 1) + (7 + 1) 



(37) 



08 + 1) + (7 + 1) 

One sees that while /12 stays constant in the drift spaces, \i\ 
and /13 change; while \i\ changes monotonously, {12 and /13 
can be locally increasing and decreasing; and etc. Besides 
that, let us note that the multiplication of any matrix T m 
from the left by a constant rotation matrix does not change 
the corresponding phase advance. For example, for an arbi- 
trary constant angle i/j the T-algorithm associated with the 
matrix 

cos(V0 sin(^) 
- sin(?/>) cos(V') 



•T, 



(38) 



will produce the same phase advance as in the original 
Courant-Snyder case, but, in general, the triangular form 
of the transition to the normalized variables ((26), which 
is also important for ID theory, will be lost. So, it seems 
that the correct way to extend the Courant-Snyder choice of 
the ID phase advance to the multidimensional case without 
losing what else important comes through the direct gener- 
alization of the ID T-algorithm defined by the matrix T\ to 
many dimensions. 

The matrix T\ is a lower triangular matrix with positive 
diagonal elements and as its multidimensional analog we 
will take the lower block triangular symplectic matrix 







(39) 



where u and w are, respectively, a symmetric and an in- 
vertible n x n matrix. Note that any lower block triangular 
symplectic matrix can be represented in this form with the 
proper choice of the matrices u and w. 

Proposition 2 The beam matrix E can be diagonalized by 
congruence using symplectic block triangular transforma- 
tion of the form ( 1391 ) if and only if 

{^J) 2 qp — ^qq ^qp ~ ^qp ^qq = Oj (40) 
and the condition ( 1401 ) is invariant under linear symplectic 
transport (0 of the beam matrix E if and only if the matrix 
E is the matrix with equal eigenemittances. 

This proposition tells us that though not only E matrices 
with equal eigenemittances can be diagonalized by the ma- 
trix d39l , the consistent theory of diagonalization by the 
lower block triangular symplectic matrices can be created 
only for the beams with equal eigenemittances. 

Substituting representation (|39l into relation (IT4t we ob- 
tain the following equations for the determination of the 
matrices w and u as functions of the given matrix Twiss 
parameters P, a and 7 

P = ww T , (41) 

a — w u w^ 1 , (42) 

7 = w~ T (I + u 2 ) w- 1 . (43) 

The general solution of the equation (|4TT > can be written as 
w = wr, where w is any particular solution of this equa- 
tion and r is an arbitrary n x n orthogonal matrix. When 
some solution of the equation (|4TT > is chosen, then the ma- 
trix u is uniquely defined from the equation (l42l and is 
symmetric, and the equation d43l is satisfied automatically 
due to relation ( TT6b . Thus the only remaining uncertainty in 
the choice of the T-algorithm lies in the nonuniques of the 
solution of equation ( |4TT >. Unfortunately, this uncertainty 
cannot be resolved simply by the requirement to recover 
the Courant-Snyder choice in the limit of multidimensional 
uncoupled beams. For example, both, w taken as unique 
Cholesky factor of the matrix p and w = P 1 ^ 2 , where p 1 / 2 
denotes an unique positive definite symmetric square root 
of the matrix P, will satisfy this requirement. But if we will 
take into account that the Courant-Snyder matrix T\ has 
positive diagonal elements and will reformulate this prop- 
erty in the form that lxl diagonal submatrices of the ma- 
trix Ti are symmetric positive definite matrices, then the 
choice for w becomes unique. One takes w — p 1 ! 2 and 
obtains for the matrix T the final (and familiar from ID 
theory) form 

jS-Va 
p-^a P 1 ' 2 



T = 



(44) 
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